The Excitation of Planar Dielectric 
Waveguides at p-n Junctions, I 

By J. McKENNA 

(Manuscript received April 26. 1967) 

The fields excited within a planar dielectric waveguide by an externally 
incident electromagnetic field are studied in this paper. The dielectric 
waveguide Jills the half space z > 0, while the halj space z < is air. 
The waveguide is formed by a nonuniform, anisotropic, nonabsorbing, 
dielectric medium. Different choices of the dielectric tensor for this medium 
yield different waveguides. Certain models which are particularly relevant 
to electro-optic diode waveguides and laser diode amplifiers are studied 
in some detail. An arbitary incident field will, in general, excite not only 
a finite number of propagating modes, but also a background of continuum 
modes. Integral representations of the total transmitted field within the 
waveguide as well as of the reflected field are obtained. The representation 
of the total transmitted field can be decomposed into a finite sum of discrete 
propagating modes, a continuum propagating field, and an evanescent 
field. Explicit evaluation of the fields depends on the solution of a pair 
of integral equations. In practice, the dielectric tensor of the waveguide 
differs but little from the dielectric constant of the surrounding material. 
An approximate solution is found for this case, and numerical results 
will appear in a following paper. 

I. INTRODUCTION 

Recently there has been great interest in the guiding of light by the 
p-n junction region in certain piezoelectric semiconductors, for it has 
been noted that the Pockels effect due to the electric field within the 
p-n junction can be used to modulate light which propagates parallel 
to the junction plane. 1 ' 4 This effect was first observed, and has been 
most intensively studied, with visible light in GaP junctions, 1 but it 
has also been observed with infrared light in GaAs junctions. 1,4 

All treatments of the effect so far have assumed that the p-n junc- 
tion region, which has a higher dielectric constant than the surround- 
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ing, normal GaP, behaves like a dielectric waveguide. 1-8 A detailed 
analysis of this waveguide would require a knowledge of the optical 
properties in the neighborhood of the junction. However, since these 
properties change significantly in a fraction of a wavelength, it is ex- 
tremely difficult to investigate them individually by experimental 
means. In order to get around this difficulty it has been necessary to 
adopt an indirect approach based on analyzing a number of different 
mathematical models and comparing their predictions with experiment. 

As part of this program Nelson and McKcnna have investigated 
the possible discrete modes which can propagate in a number of dif- 
ferent models and have studied in considerable detail the properties of 
the lowest-order mode of each polarization. Recent experimental work 
has made it increasingly clear, however, that a knowledge of the dis- 
crete modes alone is not enough to provide an understanding of these 
p-n junction dielectric waveguides. This is because a beam of light, 
when focused on the face of a junction waveguide, excites within the 
waveguide not only a finite number of discrete modes, but also a back- 
ground of continuum modes. In many cases this background light is 
intense enough to mask important features of the discrete propagat- 
ing modes. Thus, unless an understanding of this background light is 
available, the task of comparing the predictions of different mathe- 
matical models with experiment is almost impossible. An understand- 
ing of the electromagnetic boundary value problem involved also has 
great relevance to understanding what happens when light is intro- 
duced into a laser diode amplifier. 

The purpose of this paper is to study in some detail a class of math- 
ematical models of the excitation of dielectric waveguides. These mod- 
els are simple enough so that the mathematical analysis can be per- 
formed and the background light can be investigated carefully. At the 
same time, it is felt that the models are realistic enough so that their 
predictions can be compared with experiment. 

The models can be described as follows. The waveguide consists of 
the half space z > 0, as shown in Fig. 1, while the region z < is air. 
The waveguide itself is assumed to be formed by a nonuniform, aniso- 
tropic, nonabsorbing dielectric. The components of the dielectric ten- 
sor are functions of the coordinate x only, and for each value of x the 
dielectric tensor is diagonal in the fixed coordinate system shown in 
Fig. 1. As an example, for the GaP electro-optic diode modulator stud- 
ied in NM this corresponds to the cases where the junction field is in 
the [111] or [100] directions. Each such model is determined by its 
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Fig. 1 — Symmetric step model illustrating the coordinate system used in all 
the models. The dielectric tensor is always diagonal in this fixed coordinate sys- 
tem. 



dielectric tensor whose diagonal elements we will denote by K„(x) 
(n = x,y,z). 

It was shown in NM that the amount of absorption encountered in 
GaP electro-optic diode modulators was too small to affect significantly 
the shape of the modes. It is, therefore, felt that the study of absorp- 
tionless models here is well justified. It was also shown in XM that 
the detailed analytical form of the functions K„(x) is not important 
when only the lowest-order discrete mode of each polarization can 
propagate. The most important features of the discrete modes can be 
determined by studying models for which the functions K„{x) are 
step functions (piece-wise constant). Although it is possible to carry 
out a good deal of the analysis without specifying the functions K„(x), 
the final detailed results naturally depend on the choice of K„(x). We 
shall concentrate here on two models, the symmetric step model and 
the asymmetric step model. The symmetric step model is defined by 
the equations 6 

/v„,(.r) = K m , | x | < w (1) 

= K a , | x | > w (2) 

and the asymmetric step model is defined by the equations 

KJ,r) = K m , | x | < w (3) 
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= A', , x < -w (4) 

= K a , x > w, (5) 

where K, < K lt and K m > Kj ^ 1, m = x, y, z, j = 0, 1, 2 (see Fig. 

2). In the case of the GaP electro-optic diode modulators there are 

relations of the form 



K m = n 2 (l + 8,„), (m = x, y,z) 



(6) 



A' = /r(l - A), K i = n 2 (l - A,), (j = 1, 2). (7) 

In (6) and (7) n is the index of refraction of normal GaP, the quan- 
tities 8 m are linear in the junction field (the linear electro-optic effect), 
andO ^ |8,„| < A « 1. 

In Section II we will write down general integral representations for 
incident waves in the region z < 0, as well as integral representations 
for the resulting reflected and transmitted fields. These, integral rep- 
resentations will involve a number of unknown functions. Some of 
these functions are determined directly from the structure of the wave- 
guide and are independent of the incident field and the boundary con- 
dition at 2 = 0. The remaining unknown functions are determined by 
the incident field and the boundary conditions at z = 0. We show that 
these functions satisfy a set of linear integral equations. The results 
of Section II are independent of the specific form of K„,{x) and the 
incident field. In Section III we explicitly calculate the unknown func- 
tions which depend only on the structure of the waveguide for the 
symmetric and asymmetric step models. In Section IV we obtain ap- 
proximate solutions of the integral equations for a special class of 
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Fig. 2— (a) The function K,„(x) for the symmetric step model, (b) The 
function Km(x) for the asymmetric step model, 
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waveguide models. The remaining unknown functions are determined 
for these models in terms of the incident field. In a second paper on 
this subject we will give asymptotic expansions and numerical results 
for the fields within the waveguide for the symmetric and asymmetric 
models when they are excited by a Gaussian incident wave. 

II. A GENERAL DESCRIPTION OF THE FIELDS 

In this section we study formal solutions of Maxwell's equations 
which describe an incident wave in the region z < moving to the 
right and striking the waveguide from the left, a reflected wave in the 
region z < 0, and a transmitted wave in the region z > 0. The fields 
are assumed to be monochromatic and independent of the coordinate 
y. We write for the total electric and magnetic field vectors 

E(.r, z, I) = Re (e(.r, z)e iul ), H(.r, z, t) = Re (h(.r, z)e iu "), (8) 
and for the total electric displacement and magnetic induction vectors 

D(.r, z, t) = Re (d(.r, z)e iul ), B(x, z, t) = Re (b(.r, z)e^'), (9) 

where Re denotes the real part and u = 2tt/ is the angular frequency 
of the radiation. Then Maxwell's equations are 

V X e = -tub. V-d =0, 

V X h = j'«d, V-b =0. 

From our assumptions about the model, the constitutive equations can 
be written as 

b = M (,h, d = e K e, (11) 

where e and Mo are, respectively, the permittivity and permeability 
of free space. The dielectric matrix K = K(.r, z) is the unit matrix 
for z < 0, and for z > it is a diagonal matrix whose diagonal elements, 
K n (x) in = x, y, z), are functions of x only. It is a straightforward 
matter to show that any solution of Maxwell's equations satisfying 
the above assumptions can be written as the linear combination of a 
TE solution and a T M solution. We consider these solutions separately. 

2.1 TE Fields 

Wo first look for TE solutions having the form 

e(.r, z) = [0, e„(x, z), 0], h(x, z) = [h r (x, 2). 0, h : (x, z)}. (12) 
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In the region z < 0, e y must satisfy the Helmholtz equation 



^1 + ^5 + k\ = 0, (13) 

dx dz 

where the free-space wavenumber k is defined by 

k = co(eoMo) 1 = 2r/X 

and A is the free-space wavelength. The total field in z < is the sum 
of the incident field ej° and the reflected field ej r) and both ej° and c l u T) 
are solutions of (13). In the region z > there is only the transmitted 
field which satisfies the equation 

ft + ft + k 2 K v (x)e u = 0. (14) 

A solution of (13), which can be found by separation of variables, 
and which describes a general incident field due to sources in z < 
at a finite distance from the plane z = 0, is 

el'\x, z) - gj f SJ°(Z) exp { -*Q0)« ~ ^) «i < 15 ) 

where 

m = + VF^-~T, 1 1| ^ * (16) 

= -iVF^¥, I / 1 ^ ft- 

The components of the magnetic field vector can be obtained with the 
aid of Maxwell's equations by differentiating (15). Let 2(z ) denote 
the strip — oo < x < oo, ^ ?/ = 1, ly in S in tnc plane z = 2 »- Tnen 
the time averaged power incident on 2(z), * = 0, is independent of z 
and is 

p, = -i Re f e'°(.r,z)/i< (:r.,2)* tf.c 

J-oo 

= (4™*,)"' f Vfc^T I 8{ o (0 | a A, (17) 

where * denotes complex conjugation. We will assume that 

P I 6j o (0 I 1 cU < °° and f° | fl(/) I I fij°(0 P « < °° • 

(15) is to describe an incident field due to sources at z = - «*, then 
it is easy to see that we must have fij n (0 = 0, 1 1 1 > ft. Since the incident 
field must be specified, it will always be assumed that & l y l) (l) is known. 
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A solution of (13) describing a general wave reflected from the wave- 
guide surface z = is 



[ r) (x, z) = — j sj r, (0 exp \iil(l)z - ih\ dl. 



(18) 



We will always assume that the source of the incident radiation is 
perfectly absorbing so that ej°(x, z) -f ej r, (.r, z) is the total field in 
the region between the source and the surface of the waveguide at 
z = 0. It will be seen that because of the boundary conditions at z = 0, 
S^ (I) generally does not vanish outside some finite I interval. Because 
of the factor exp [iQ(l)z\, that part of the integral in (18) between 
the limits -/.■ and /,-, j k _ k { j dl, represents a traveling field, while 
the remainder of the integral represents an evanescent field which 
damps out very rapidly in the negative z direction. The time averaged 
power reflected back through the strip 2(2), z ^ 0, is 

P r = (4™^)" f Vk* - I 2 I 6* r) (i) | 2 dl. (19) 

We now turn to the transmitted field. We use the method of separa- 
tion of variables, and we seek transmitted waves which are linear 
superpositions of solutions of (14) of the form 



.<n, 



1 (.r , 2) % e y (x) exp | - i V - vz } . (20) 

In (20) p is a real separation parameter, and if v > 0, \/— ~v = —i\/p. 
If (20) is substituted into (14) we get the eigenvalue equation 

—^ + (A-7v„(.r) + v)e a = 0. (21) 

Equation (21) defines a singular, self-adjoint, second-order boun- 
dary value problem on the interval — 00 < x < 00. The theory of this 
equation is well known, and we refer the reader to Coddington and 
Levinson 7 for a detailed treatment. We give a summary here of those 
properties of such equations which we will need. 

For all the models under consideration, the functions K m (x) are 
positive, bounded functions, which are bounded away from zero, and 
which are differentiable except for at most a finite number of step dis- 
continuities. Equation (21), therefore, defines a problem which is called 
limit-point type at both plus and minus infinity. This means that for 
arbitrary, complex v, (21 ) possesses exactly one solution (up to a con- 
stant factor) which is square integrable over < x < 00, and exactly 
one solution which is square integrable over — 00 < x < 0. 
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For a given real number u, let <p,(.t, v) and <p 2 (x, v) be the two solutions 
of (21) which satisfy the conditions that <p,(x, v) and ^{x, v) be con- 
tinuous and which satisfy the initial conditions 

*((>,*) -1, ^(0,,)=0, (22) 

ft(0,0-0, ^(0,^ = 1, (23) 

where ' = d/dx. Equation (21) also determines a 2 X 2 matrix-valued 
function p(v), — oo < v < w, having the following properties: (i) piy) 
is Hermitian { Pik {v) = p£ (<>)). (*0 p(") _ p( w ) is P ositive semidefinite 
if j/ > m- (w) P,t(") is of bounded variation on every finite interval. 
The matrix p{v) is called the spectral density matrix and its construction 
is outlined in Section III. Then if j(x) is any square integrable function 
(/-•I /(*) I*** < °°)j we can define two transforms of j(x), g t {v) (;'=1. 2), 
such that 



lim f Z L(v) - P /(*)*>,(*, r) cfc 

i,-oo •'-oo ;.*-! >. ■'-Jt 

' W v) "* /" ffrWfc "> **/ *"*W - 0- (24a) 

'This is referred to as convergence in the mean with respect to the 
measure p(v), and in the manner of Fourier transforms of £ 2 functions, 
we write 

g,{v) = /" fGbitot ") <** (i = 1. 2). (24b) 

In terms of these transforms, the Parseval equality 

P|/(x)| f <fc- E( 9i{v)*g k {v) dp ik (v), (25) 

J-M |.*»1 "-80 

and the expansion 

/(*) = E f " ftfc. r)ff*00 dp,*(") (26) 

/.* = 1 J-x 

are valid. Equation (26) is defined in terms of convergence in the 
mean. The set of real points v at which the functions Pjk (v) are noncon- 
stant is the spectrum of (21). The set of points where any PJ M is 
discontinuous is the point spectrum and for each such value of v, (21) 
has exactly one square integrable solution. The continuous spectrum 
is the set of points of continuity of P {v) which are in the spectrum. In 
Section III we will exhibit the spectral density matrices for two im- 
portant models. 
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We can now write down a formal expression for the transmitted 
field: 

(•l n (.v,z)= £ f cxp\-iV-uzW l (x.p)g,(p)d Pik (v). (27) 

l.t-l >'-x 

The two initial value solutions <p f (x, v) (j = 1, 2), as well as the func- 
tions pjkiy) (j, k — 1,2) are determined, independently of the bound- 
ary conditions at 2 = 0, by (21) and we can assume that they are 
known. The two unknown functions gj(v) (j = 1, 2) in (27) are de- 
termined by the field at 2 = 0, since with the aid of (24) we can write 

0/00 = [ <%°(.r, 0) Vi (x, v) dx. (28) 

It is clear that because of the factor exp \—i\/—vz\, the parts of the 
integrals /"„ in (27) represent the propagating portion of the trans- 
mitted field, while the parts /" reprent the evanescent portion of the 
transmitted field. With the aid of the Parseval relation, (25), we can 
write down an expression for the time averaged power transmitted 
across any 2(2), z ^ 0, 

p, = (2^;r 2 [ v^iW*0*w<*p**0. (29) 

; <- = 1 J - f*, 



;.a-=i •>- 



We now make use of the conditions that e y (x, z) and h x (x, z) must 
be continuous at z = in order to write down a set of integral equations 
which determines &f, rl (l), <7i("), and g%(v). 

:l f &°(/) + C'aw" di = t, [° *(*, M dpM, (30) 

- » J - X I . A - - 1 •' - X 

.f f Q(0[Si°(Z) - 6i r) (0]c-" x dl 

Sir J- x 

= E ( V-^ i (x,v)q k (v)dpM. (31) 

i . k - 1 •' - x 

Although there appear to be only two equations in three unknown 
functions, because of (24) and (26). (30) and (31) are sufficient to 
determine the unknown functions. We indicate formally why this is 
true, although it will be clear from the results of Section IV that this 
scheme must be modified in specific cases. We do not go into these 
modifications, because in Section IV we use a different scheme to get 
approximate solutions. With the aid of (24b), solve (30) and (31) for 



(32) 
(33) 
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0;(v), giving the four equations 

g fy) = r nb, v) dx ± r [tf } (o + c(o]<r*" di, 

.'_„ -sir .'-oo 

V^wgfy) - f <Pi(*, p) dx ± f OC0[8i o (Q - 8i r, (0>"' fa «i 

J_oo — T J -to 

for j = 1, 2. On eliminating g t fy) between the.se equations we get the 
two equations in the unknown 8* r> (/)- 

f ^(x, u) dx ~ f (V=J + Q(0)8i r) (0e-"* <//, 

= f ^(x, r) rf.c ~ P (- V-T + OC0)8i°(0«" <ta «, 

for ; = 1, 2. Now from (26) we can write f{x) = h(x) + / 2 (*) where 

Ux) = t, P *(*, r)fcM <*P,*to (fc - 1, 2), (35) 

1=1 •'-00 

P fJMnto, v) dx = 8 n .gM (j, k - 1, 2), (36) 

and 5,-fc is the Kronecker delta function. It is this decomposition of 
an arbitrary /(.r) into components lying in the two subspaces spanned 
by tpi(x, v) and <p 2 (x, v) which is reflected in the two integral equations 
(34). The solution of (34) with given j yields the component of the 
reflected field lying in the subspace spanned by the corresponding 
tpiix, v). Let Bg'CO, 3 = h 2, denote the two solutions. Then 8j r> (Z) = 
& ( „i(l) + 8S } (0 describes the total reflected field. With this result 
g .{y) (j = \ t 2) can be obtained from either (32) or (33). We have been 
unable to obtain exact solutions for the integral equations (30)-(31) 
for any of the models considered here. However, in Section IV approx- 
imate solutions are obtained for certain situations of interest, 

2.2 TM Fields 

We next seek TM solutions of Maxwell's equations of the form 

e(x, z) = (e x (x, z), 0, e,(x, z)), h(x, z) = (0, h,(x, z), 0). (37) 
In the region z < 0, h„ must satisfy (13) . In the region z > 0, h y must 
satisfy the equation 

£ {(l/K.M) f} + l (WW f} + W. - 0. (38) 
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Just as for the TE fields, a general incident field due to sources in 
z < at a finite distance from the plane z = is 

&"'(£,*) = 5- f 3Cj°(Z) exp {-t'fi(Z) 2 - iTs) dZ. (39) 

The time averaged power due to this wave which is incident on 2(z), 
z ^ 0, is 

Pi = J Re r ei°(3, •)*i Cr l 2)* dx = (4™,,)-* [ fi(Z) | 3Ci°(Z) | 2 dZ. 

•'-CO J-k 

(40) 

We assume that /?„ | 3Cj°(0 | 2 dZ < co and /-„, | fi(Z) | | acj°(/) | 2 dZ < 
As for the 77? field if the sources of the 7W field are at z = - 00 
then 3Cj°(Z) = 0, | Z | > fc. Furthermore, it will always be assumed 
that acJ°(Z) is known. 
A solution of (13) describing a general reflected wave is 

hi'\x, Z) = ~ [ 3C l u r) (l) exp {iQ(Z> - t'Zz} dZ. (41) 

Just as in the case of the TE field, h v T) (x, z) can be split into a prop- 
agating field and an evanescent field. The time averaged power re- 
flected back through the strip Z(z), z ^ 0, is 

P, = (4fl Weo )- 1 J h fi(Z) I 3C: r) (Z) | 2 dZ. (42) 

The transmitted field is again treated by separation of variables, 
and we write 

h u l) (x, z) « h y (x) exp { -i'V-^2] . 
Then /*„(.!•) satisfies the eigenvalue equation 

KX-c) ^ |(l/A',(.r)) *j| + (*■*,(?) + ,)A. = 0. (43) 

Equation (43) is not in the canonical form of a self-adjoint boundary 
value problem. However, if we make the change of variables 

u = f [K x (£)\~ l dl, (44) 

(43) is transformed to the equation 

£ [{KMK.(u)\-' fj] + PKM + v)h u = 0. (45) 
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This equation defines a self-adjoint boundary value problem, 8 and 
even though the function [K x (u)K,(u) } _1 may have step discontinuities, 
the techniques of Ref. 7 can be shown to be still valid. Equation (45) 
is limit-point at u = ± » , and so on transforming back to the variable 
x, the following statements can be made. 

For a given real number v, let ^,(x, v) and f 2 (x, v) be the two solu- 
tions of (43) which satisfy the requirements that 

*,(?,*) and {K t {x))- l ^{x,v) 
be continuous for all x, and which satisfy the initial conditions 

*i(0, v) = 1, (l/£.(0))#(0, ,) = 0, (46) 

fc(0, ") - 1, (1/tf .(0))#(0, v) = 1. (47) 

Equation (43) determines a 2 X 2 spectral density matrix a(v) whose 

construction is given in Section III. If f(x) is any square integrable 

function of x, we define two transforms of f(x) , 

*f« - P /(*)*.(*, ")fMr' ** o* - i. 2)- < 48 ) 

•'-00 

where equality in (48) is defined in terms of convergence in the mean 
with respect to the measure o-(v). In terms of these transforms, the 
Parseval equality 

f \j(x) | 2 [KM)' 1 dx= E f MtfW <M>), ( 49 ) 
and the expansion 

/(*) - tr^W)^ (50) 

are valid. The last equality is again defined in the sense of convergence 

in the mean. 
We can write down a formal expression for the transmitted field 



2 /«eo 

hy\x t z) - 2 / ex P {— *V— y*iiM* f v)h k (v)d<r ik (p). 

i.k = \ •'-oo 



(51) 



The two initial value solutions fyix, v) (j = 1, 2), as well as the func- 
tions a jk (v) (j, k = 1,2) are determined, independently of the bound- 
ary conditions at z = 0, by (43) and we can assume that they are 
known. The two unknown functions hj(v) (j = 1, 2) in (51) are de- 
termined by the field at z = since, with the aid of (48) we can write 
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M0 - f K n (x,0)U*,v)iK z (x)}- 1 dx. (52) 

J-00 

With the aid of the Parseval relation, (49) , we can write down an ex- 
pression for the time averaged power transmitted across any 2(2), 
2^0: 

P, - (2W60)" 1 £ f V^hM*h k (u)dcr ik ( v ). (53) 

We can now make use of the conditions that e x (x, z) and h„(x, z) 
must be continuous at z = in order to write down a set of integral 
equations from which 3C< r) (0, h v (v), and h 2 (v) can be determined. 

£ j [«i°(0 + ^"(QK "* dl= t T Ux, v)h k {v) daM, (54) 

~£ QiDl^d) - 3cl r \l)] e - ilz dl 



2 .,» 

j.*-i J-00 



(55) 



Just as in the case of the TE field, the solution of (54) and (55) re- 
duces to the solution of the two integral equations 

f *,(*, »)dx~ f j V^/K,(x) + n(l)M r) (l)e- ilx dl 

J — CO £" J — x 

= f *,(», v) dx ~- f j - V~v/K x {x) + Q(Z) } 

.ac{°(0«-' ta ctt, (j - 1, 2). (56) 

III. THE SPECTRAL DENSITY MATRIX FOR SEVERAL MODELS 

3.1 General Outline of the Construction 

In Section II it was shown that the determination of the transmitted 
field for a given model depended on a knowledge of the initial value 
solutions <p,(x, v) and #,(x, v) (j = 1, 2) and the spectral density ma- 
trices p{y) and a{v). In this section we study these functions in some 
detail for two simple but important models, the symmetric step model 
and the asymmetric step model. These calculations illustrate the 
technique for treating the whole class of piecewise constant models. 

We first outline the general construction of the spectral density 
matrices. 7 The solutions of (21) have the property that the functions 
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<Pi(x, v), (p'jix, v) (j = 1, 2) are entire functions of v for each fixed x, 
when v is a complex variable. The first step is to determine the two 
functions of v, m„0>) and m_„(iO such that when Im v > 0, ^>i(x, v) + 
m m (i>)(p 2 (x, v) is a square integrable function of x over [0, °°] and 
<Pi(x, v) + m-a,(v)<p 2 (x, v) is square integrable over [-<», 0]. The ele- 
ments of the spectral density matrix are then given by the formula 

1 f ' 

Pikiy) - pM = lim - / Im M ik (n + ie) dv (57) 

where p. and v are real, Im denotes the imaginary part, and for arbi- 
trary complex v 

MM = (m-GO - mjy))~\ ( 58 ) 

Af 19 (v) = M 21 (.) = \(m-Jy) + m 00 W)(m_ 00 W - mJy)T\ (59) 

M 22 (.) = m.»wi.W(wi_» - mJy)Y l . (60) 

Equation (57) defines p ik (y) uniquely at points of continuity up to 
an arbitrary, additive constant. The functions M ik {v) (j, k = 1, 2) 
are meromorphic if Im v j* and all their real poles are simple. The 
point spectrum consists exactly of the points which are real poles 
of one of the M ik {v). There are at most a countable number of such 
points. Let v Q be a real pole of M ik (v) and let a ik be the residue there, 

MM=-^-+ ••• • (61) 

v — v 

Then it follows from (57) and (61) that 

P,t(Po + 0) - Pik (v - 0) - -Re (a it ). (62) 

If y is not a pole of any M ik {v), and Im M ik (u ) ^ for some (j, k), 
then 1*0 is a point of the continuous spectrum and 

d Pi M =-ImI,W. (63) 

IT 

If v is not a pole of any M jk (v) and Im MjM = for all (;*, k) in 
some neighborhood of v , then v is not in the spectrum and 

dpM = (j, k = 1, 2) (64) 

in a neighborhood of vo. 

3.2 TE Fields for Symmetric Step Model 

We now apply these formulas to the symmetric step model for the 
case of the TE field. The functions K n (x) (n = x, y, z) are defined by 
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(1 and 2). Equation (21) has constant coefficients in the two regions 
\x\ < w and \x\ > w. Since e„(.r, z) and h,(x, z) must be continuous 
at x = ±w, the desired solution of (21) must be continuous and have 
a continuous derivative. We have 

¥>i(z, v) = COS (co u x), | x | ^ w (65) 

= cos (u u w) cos (« (| x I — w)\ 

- k/w ) sin (a),™) sin jo> (| *|-w)j, | a: | ^ w (66) 

¥*G&, v) = (l/w B ) sin (o>„a;), \ x \ ^ w (67) 

= (l/a>„) sin (a) v w) cos {w (x — w) } 

+ (l/co„) cos (c^w) sin \a) (x - w)}, x ^ w (68) 

^(z, J') = -<f 2 (-x, v), x ^ -w (69) 

where 

a>„ = [v + fcXl 1 (n = 0, x, y) . (70) 

In (70) <o„ is defined as a single-valued function of v in the complex 
plane cut along the real axis from -k-K n to 00. That branch is chosen 
which is positive real on the upper side of the cut. Simple calculations 
now yield 

mjy) = —m-Jy) 

= {«, sin (u y w) + iu cos («»«;)} {cos (« y w) - i(fa /<a y ) sin («,tp)} -1 . 

(71) 
Therefore, 

M u (p) = -l/[4M as ( v )\ = l/{2m_»}, (72) 

il/wW = M.,00 = 0. (73) 

In order to determine the spectrum, we begin by decomposing the 
whole real axis into the union of three intervals 

h - [-oo, -k?K,-\, U = i.-k 2 K, , -k 2 K ), I 3 = [-k 2 K , 00]. (74) 

From (57) and (73) it is clear that pi2 ( v ) and P21 (v) are constant for 
all v, hence 

dp i2 (v) = dp 2l {v) =0, -co^j/^co. (75) 

It is easily seen that M n (v) and M-^iv) are real and have no poles or 
zeros in I x . Therefore, I x contains no points of the spectrum, and 

/>,,(") = P,-,-(-«>), dp u (y) = ve J, (j = 1,2). (76) 



1506 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1967 

In the interval I 2 , Mu(v) and M 22 {v) can each have a finite num- 
ber of poles, and from (72) it follows that the poles Mu(v) are the 
zeros of M 22 (v) and vice versa. The real poles of InW are the real 
solutions of 

<a„ sin (o) u w) + iw cos (a) u w) = (77) 

and the real poles of M 22 are the real solutions of 

cos (u)„w) - t(co /w„) sin (u v w) = 0. (78) 

For v e I 2 , *»„ is real while o is purely imaginary. If we let 

&(„) = Uy (y) } pfy) = -iu (v) = {-v - k 2 K ) h , (79) 

then (77) in the single unknown v can be replaced by the set of three 
equations 

-v = k 2 K + p 2 , -v = fc% - 6 2 , b tan bw = p, (80) 
in the two positive real unknowns b and p and the original unknown 
v. Similarly, (78) can be replaced by the set of equations 

-v = k 2 K + p 2 , -v = k 2 K„ - b 2 , b cot bw = -p. (81) 
These equations are well known and their solutions have been deter- 
mined. 6 - 9 The set of equations (80) has a finite number of real solu- 
tions and always has at least one solution for all positive values of the 
parameters, w, k, K y - K . These are the even modes of NM. We 
denote corresponding values of v by viy, j = 1, 2, • • • , Ri. The set of 
equations (81) also has a most finite number of solutions, although if 
{wk) z X (K y — K ) is small enough it has no real solutions. These are 
the odd modes of NM. We denote the values of v corresponding to these 
roots v 2h j = 1, 2, ■ • •, R 2 . The points v lh v 2h which are all in the in- 
terval I 2 , comprise the point spectrum of (21) . Let 

8p(v) = lim \ P ( V + e) - p(v - «=)}. (82) 

«-»+o 

Then with the aid of (62) it is easy to show that 

SpnKO = p(vu)/{l + wpbu)), WO = 0, i = 1,2,--,B X , (83) 

Spn("2,) = 0, W"2,) = b 2 (v 2j )pM/{l + wpM), 

j = 1,2, ••• ,R* . (84) 

With the aid of (65) through (69) and (77) through (79) it is readily 
shown that 
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9i(x, t>ij) = cos (b(v u )x), I x I ^ w (85) 

= cos (b(*,f)iiO exp {p(vi,)(«J - [ .r |)|, | x | ^ w (86) 

f2(x, v 2i ) = \l/b(v 2i )} sin (b(v ai )x), \ x \ £ w (87) 

= (l/6(^ 2 ,)} sin (b(p 2j )w) exp \p(v ti )(w - x)\, x £ w (88) 

<Pa(z, "2,0 = — <p 2 (-x, v 2i ). x ^ —w (89) 

It is also true that 

/_«, Vi(x ' v,k)% dx = VWftOi fe - 1, 2, • ■ • , R, , j = 1,2. (90) 

The remaining points in I 2 are not in the spectrum. 

Finally, in the interval I 3 it is readily shown that M n (v) and M 2 z(v) 
have no poles. It is shown easily then that the whole interval 7 3 is in 
the continuous spectrum, and in this interval 

d P ,,{v) = p'n{v)dv O - =1,2), (91) 

where 

PiiW - IT [«*8irf (u y w) + w'o COS 2 (aj w w)]~'aj„ , (92) 

LIT 

AuOO = ^~ [w» cos 2 («„«>) + coo sin 2 (a>„w)]~'a> 2 ci;n . (93) 

In summary, the spectrum of (21) consists of the points vj k , k = 
1, 2, • • •, Rj, j = 1, 2, and the interval J s . Equation (27) for the trans- 
mitted field can be written as 

2 Ri 

ei°(z,«) = 2 £ faa(?ik)9i(vik)<Pi(x, V{k) exp {— «V— p,**] 

j • - 1 * - 1 



dp 



2 »0 

+ Z) / exp {-tV^aJ^foiOfrCiOp/'yCiO 

2 »M 

+ E / exp \-Vv'z\<p i (x, V )g i {v) P ' ii (v)di>. (94) 

»-l •'0 

The terms in the first, double summation in (94) are just the possible 
TE modes which can be excited in the waveguide. The terms in the 
second summation represent the propagating continuum field while 
the terms in the last summation represent the evanescent part of the 
transmitted field. A useful interpretation of the propagating continuum 
field can be obtained as follows. Consider within the waveguide in the 



1508 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1967 

region x < —w an incident plane wave of the form 

el 0) (x, z, v) = exp { -iV^v z - iu (v)x) , (95) 

so that if 6 is the direction of propagation of this wave (measured 
clockwise from the positive z axis) , then 

cos 6 = V^v/kVWo , sine = u (v)/kVK~ . (95) 

On striking the region of higher dielectric constant, \x\ < w, part of 
this wave will be reflected and part of it will be transmitted through 
the region \x\ < w. Denote by x + ( x > z > ") this total electromagnetic 
field set up by the incident wave, (95) . Similarly, denote by X - ( x , z > *) 
the total electromagnetic field set up by the incident wave in the region 
x > w 



e { v °\x, z, v) = exp {—iV^-vz + iu Q (v)x} 



(97) 



In Fig. (3) we give a schematic description of x+ and x- • Then it can 
be shown that for —k 2 K ^ v ^ 0, 

exp {-iV^vz}(pi(x, p) = aj(v)x+(x, z, v) + &,(")x-(z, z, ") (i = !» 2 )- 

(98) 

For the above values of v the directions of propagation of the incident 
waves for x+ and X - fill the interval -t/2 ^ 6 ^ tt/2. Thus, the prop- 
agating continuum field is just a wave packet of plane waves appropriate 
to the medium defined by the dielectric tensor K„(x). 

Similarly, the evanescent part of the field can be interpreted as a 
superposition of waves bound to the surface z = and propagating in 





(b)x_ 



Fig. 3 — A schematic diagram of the plane waves appropriate to the dielectic 
medium in the symmetric step model. The wave x* is incident on the junction 
region from the positive x direction, while x- is incident from the negative x 
direction. 
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the positive and negative x directions. The distinction between the 
propagating and evanescent parts of the transmitted field is further 
shown in the expression for the time averaged transmitted power, (29) , 
which for the symmetric step model is 



Ri 

1 



P t = (20W' 1 Z £ V^ k I <7,M | 2 hiM 



+ (2^0)-' z r v^v i gM i 2 P ,» dp. 

J-l J-k'K 



(99) 



As this expression shows, the evanescent part of the field transmits no 
energy on the average. 

3.3 TM Fields For Symmetric Step Model 

The TM fields of the symmetric step model can be treated similarly. 
Equation (43) has constant coefficients in the two regions \x\ < w 
and | x | > w. Since e,(x, z) and h„(x, z) must be continuous at x = ±w, 
the solutions of (43) must be such that ^,(.r, v) and \\/K t (x)W t {x, v) 
(j = 1, 2) are continuous. We have 

fiix, v) = cos (K r u z x), | x | ^ w (100) 

= cos (K T co x w) cos {a> (| x | — w)\ 

- K«x/Q/(w 7Q} sin (K T co z w) sin {o> (| x \ - to)}, | x | ^ w (101) 

Mx, ") = {K B M sin (K r a),x), \ x | ^ w (102) 

= {/^/coj sin (K r u) T w) cos {co (x — w)} 

+ {Kq/wo} cos(K r a) x w)s'm {w (x — w)}, x ^ w (103) 

^(z, ") = -ip,(-x,v), x ^ -w (104) 

where 

K. = {K X K$, K, = {K g /K,)\ (105) 

and m w and <o are defined in (70) . Next, 

m x (v) = -m.jp) = {(a> x /K )sm (K r a>,w) + i(u /K ) cos (K r u z w)} 

■ { cos (7v r a x w) - i(K u /K co z ) sin (K r w z w) J" 1 . (106) 
Therefore, 

M u (iO = -l/!4.1/ 22 (v)} = l/!2m_ K (,)J, (107) 

JI/ l2 (iO - M 2l (p) = 0, (108) 
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and from (57) and (108) we have 

d*M = da 2l (y) =0, - * < v < oo . (109) 

The spectrum in the case of TM fields is determined in the same 
way as in the case of TE fields, and we merely state the results. There 
are no points of the spectrum in the interval h = [— oo , —k 2 K x ] , 

*M = »«(- °°). *«M = °i v £ J - G = i. 2)- ( n °) 

The interval I 2 = (-k 2 K x , -k 2 K ) contains a finite number of points 
in the point spectrum. The points of discontinuity of cm(v) are the 
real solutions of 

(u t /K g ) sin (K T u T w) + i(u> /K ) cos (K*> x w) = 0, (111) 

while the points of discontinuity of a&(v) are the real solutions of 

cos (K> x w>) - i(K>o/ff<*0 sin (Kru z w) = 0. (112) 

If we let 

bfy) = KMf), p(f) = -io»o(v) = (-p - k 2 K )\ (113) 
then (111) in the single unknown v can be replaced by the set of equa- 
tions 

-v = k 2 K n + v\ "" = k * K * ~ K * h2 / K * » bK ° tan hw = P K ' » 

(114) 

in the two positive real unknowns b and p and the original unknown 
v. In the. same way, (112) can be replaced by the set of equations 

-v = k 2 K Q + V \ -v = ***- - K z b 2 /K, , bK cot bw = -pK. . 

(115) 

The set of (114) has a finite number of real solutions and for all posi- 
tive values of the parameters K /K s , K x /K z , w, k 2 {K x - K ) there is 
always at least one solution. 6 - 9 These are the even modes of NM. The 
corresponding values of v are denoted by vxj, j = 1, 2, • ■ • , Si. The set 
of equations (115) also has at most a finite number of solutions, al- 
though if (wk) 2 (K 9 — K ) is small enough it has no real solutions. 
These are the odd modes of JVM. The corresponding values of v are 
denoted by v 2j , ; = 1, 2, • • •, S 8 . The points v v , v 2 , are the point spec- 
trum of (43) and they all lie in the interval I 2 . Furthermore, 

Struma) = S(pK)), !*>(*,) = 0, j - 1, 2, •• • , S l , (116) 
M = 0, 8<r 22 (v 2i ) = b( V2i ) 2 S(p(v 2i ))/K 2 . , j = 1, 2, • ■ • , S 2 , (117) 
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where 

m = k. p [wp + (gA -gffifig - * J '- (118) 

From (100) through (104) and (111) through (113) it follows that 
ft (a, "i.O = cos (&(v,,-).t), I a; \ £ w (119) 

= cos (b(p lj )w) exp [p(vn)(ifl — | x |)j, | a: | g£ w (120) 
Mx, v 2j ) = {K : /b{v 2l )} sin (b(p 2i )x), \x\£w (121) 

= [£,/&(i>h)} sin (b(j/ 2y ).r) exp {p(j/ 2 /)(w -a:)}, x £ w, 

(122) 



ft(z, "2,0 = -ft(-.r, v 2; ). .t ^ -w 
It is also true that 



(123) 



f fcfa ir,,) 1 ^*)}" 1 <& = l/lrnO'fi), * = 1, 2,- • -, S y , j = 1, 2. 

J-00 

(124) 

The remaining points in 7 2 are not in the spectrum. 

The continuous spectrum is the interval I 3 = [ — k 2 K , oo]. For 
points of the continuous spectrum 

dff fi (y) = o'iMdv (j = 1,2), (125) 

where 

on GO = 2~ [^nw^sin 2 (K r (O x V)) 

+ iCK\a>o cos 2 (/C r o) I «>)]" 1 A' K' I K t a; n , (126) 

0-2 2 OO = y ^« w * cos2 ( K ^ z w) 

+ #,# -to 2 , sin 2 (ff^ x uO]~Xc&>o . (127) 

To summarize these results, the spectrum consists of the points v#, 
k = 1, 2, ■ • •, Sj, ; = 1, 2 and the interval 7 3 , and the transmitted field 
can be written in the form 



h { v l \x,z) = 2 £ (rn(vik)hi(v ik )^,(x, v ik ) exp \—tV—v ik z\ 
+ Z) / exp J-iV-i'zlftte, r)Af(rMf60 dv 

j-i J-it'/ro 

2 /»ao _ 

+ Z / exp {- V> s} ft(a;, r)A,to»«W tfo. 



(128) 



1512 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1967 

Just as for the TE fields, the terms in the first, double summation in 
(128) are the possible TM modes which can be excited in the wave- 
guide. The terms in the second summation represent the propagating 
continuum field while the terms in the last summation represent the 
evanescent part of the transmitted field. Just as for the TE fields, the 
propagating part of the continuum field can be interpreted as a wave 
packet of reflected and refracted plane waves, and the evanescent part 
of the field can be interpreted in terms of surface waves at z = 0. 
Equation (53) for the transmitted energy is 



P t = (2C060)" 1 E i V^ | hi(y ik ) | 2 totM 



2 Sj 

EI 

+ (^eo)- 1 £ f v^ I *fM I 2 •foM dv - < 129 ) 

, = 1 J-k'Ko 



3.4 TE Fields For Asymmetric Step Model 

We now turn to the second of the two models which are studied in 
detail and examine the TE fields for the asymmetric step model. The 
functions K n (x) (n = x, y, z) are defined by (3) through (5). Equa- 
tion (21) has constant coefficients in the regions \x\ < w, x > w, x < 
—w, and we seek solutions which are continuous and have continuous 
first derivatives. Then 

<Pi(x, v) = cos (ayX), | x | ^ to (130) 

= cos («„to) cos \u 2 (x — w) } 

— (uyM sin (u u w) sin \u 2 (x — w)}, x ^ w (131) 
= cos (<a w to) cos |w,(.r + w)\ 

+ («,/«,) sin (uyto) sin [a^x + w)}, x ^ — w (132) 
<p 2 (x, v) = (l/co„) sin (w„u>)> | x | ^ to (133) 

= (l/«„) sin (coyiv) cos [a> 2 (x — w)\ 

+ (l/« a ) cos (w„w) sin {w a (a; - w)}, x > w (134) 
= — (1/oia) sin («„«>) cos {w,(a; + if)} 

+ (l/coi) cos (<o w w) sin { Wl (.r + «;)}, a: ^ -w, (135) 
where 

«.00 = (" + W (n = 1, 2, a, t/). (136) 
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As before <% is defined as a single-valued function of v in the complex 
plane cut along the real axis from —k 2 K n to oo. Then 

mjy) = {w v sin (a v w) + iu a cos (co u w)} 

• {cos {u v w) — i(u) 2 /o} u ) sin (io v w)}~ 1 , (137) 
m-„(v) = — {o3„ sin (a>„w) + ?a>, cos («„«>)) 

• {cos (u y w) — i(«i/«„) sin (« B w)}~ 1 . (138) 
From (58) through (60) and (137) through (138) we obtain 

M ik (y) = N ik {v)/D{v) fj, k = 1, 2), (139) 

where 

#ii(»0 = -£[(1 ~ w.coz/co') + (1 + co^/co*) cos (2u v w) 

- &{(«, + aj 2 )/w„} sin (2co v w)], (140) 
iV 12 W = N 21 (p) = (i/2)( Wl - «„), (141) 

^22^) = 2K^y _ WiOJ 2 ) — (<4 + O)^) COS (2w„U;) 

+ tWy^! + co 2 ) sin (2(0^) } , (142) 

D(v) = (u u + u^a/ai,,) sin (2w„w>) + ?'(w, + a> 2 ) cos (2^,,™). (143) 

To determine the spectrum we note first that in the interval Ii = 
[-co, -k-K y ], the functions M jk (v) (j, fc = 1, 2) are analytic and 
real. This interval, therefore, contains no points of the spectrum and 

dpM =0 (j, k = 1,2), vtl, . (144) 

The only real poles of the functions M jk { v ) are in the interval I 2 = 
(-k 2 K u , -k^Ki). These poles are the real solutions of D(v) = 0. In 
I 2 , o> u is real while w : and w 2 are purely imaginary. If we let 

b(y) - u u { v ), Vn ( v ) = -iu n (p) = (-p - k 2 K„y (n = 1,2), (145) 

then the equation D (v) =0 is equivalent to the set of four equations 

-v = k 2 K, + pi , -v = k 2 K 2 + pi , -v = k 2 K v - b 2 , (146) 

tan 2bw = [pjb + p a /b}/[l - (p 1 /b)(p 2 /b)} , 

in the three positive real unknowns b, p lf p 2 and the original unknown 
v. These equations and their solutions have also been studied in de- 
tail. 6 - 6 In order that (146) have a solution, it is necessary and sum- 
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cient that 

K u >K n (n=l,2), (147) 

2wk(K v - /£,)' > tan" 1 {(£, - K 2 )/(K y - KJ}*. 
If conditions (147) are satisfied, D(v) =0 has a finite number of real 
solutions, VJ , j = 1, 2, •• •, R which all lie in the interval I 2 . This is 
the first significant difference between the symmetric and asymmetric 
step models. The symmetric step model always has at least one point 
in its point spectrum while the asymmetric step model may have no 
point spectrum. 

We can write, assuming that (146) and (147) are satisfied. 

SpM = -NM/Wd, i. * - !• 2. 1 = 1,2, ■■■ ,R, (148) 
where D'{v) = {d/dv) D (.<) . If we make use of (145 1 , it is easy to show 
that 

[SpM)* = &Pn(vi)*P22(vi), I = 1, 2, ••• ,R. (149) 

Neither of the functions ^(x, v,) or <p a (x, p s ) is square integrable over 
- oo < x < a. for ;' = 1,2, • • • , R. However, because of (149), they 
appear in (27) for e ( u n (x, z) only in the combination 

$(x, v,) = Vsputyj) <pi(x, v,) 

+ [SpM/V&p^MMx, v,), j = 1, 2, • ■ • , R. (150) 
If we define 

*o(», vt) = VdpiM cos (b(vi)x) 

+ I ipiM/VtpxM Hvi)) sin (fi(v t )x), (151) 
then because of (146) 

*(*, v t ) = $o(x, v,), \x\£w (152) 

= $ (w, v,) exp \ V2 {v,) (w - x)\, x = w (153) 

= $ Q (-w, Vi ) exp {p,(v,)(w + x)}. x ^ -w (154) 

Thus, the functions $(rr, vj) are square integrable, and, as we shall see, 
are just the possible propagating modes in the wave guide. The re- 
maining points in the interval I 2 are not in the spectrum. 

The remainder of the real axis, the interval -1<?K X = v ^ «>, forms 
the continuous spectrum. To show this, consider first the interval 
I 3 = [-k 2 K 1} -JrK 2 ]. In h , »» and co, are real, while o> 2 is purely 
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imaginary. The functions M ik (p) have no poles in I 3 and their imaginary 
parts are not zero. We introduce the notation 

u„(y) = b(y), t»i(f) = pi(v), w 2 (»0 = ip 2 {v), v e I, . (155) 
Then we can write 

dpM = \ \Px<?)/*b)\rfy)rfr) dv {j, k = 1, 2), (156) 

where 

ri(v) = cos bw + (p 2 /b) sin bw, (157) 

r 2 {v) = b sin bw — p 2 cos bw, (158) 

A(v) = {6 sin 2bw — p 2 cos 2bw\ 2 

+ {(p,p 2 /&) sin 2&W + p x cos 2&w} 2 . (159) 

For v c 7 3 it is clear from (131), (134), and (155) that <p^(x, v) and <p 2 (x, v) 
both grow exponentially as X — » +°°. However, from (156) we see 
that in (27) for e[ l) (x, z), the functions <p,(z, v) (j = 1, 2) appear only 
in the combination 

A(.r, v) = r x {v) Vl {x t v) + r 2 {v)<p 2 {x, v) (160) 

when v e I a . However, 

A(.t, r) = cos {6(.t — w)\ — (p 2 /b) sin {&(* — w)}, \x\ £ w (161) 

= exp {p 2 (w — x) ) , a: ^ w (162) 

= (cos 26w + (p 2 /b) sin 26w>) cos \p t (x + w)} 

+ (l/p,)(6 sin 2bw — p 2 cos 2bw) 

•sin {p^x + «>)}, .r ^ -w. (163) 

Equations (161) through (163) represent the second important differ- 
ence between the symmetric and asymmetric step models. In the sym- 
metric model all the components of the continuum field are oscillatory 
functions of x on both sides of the waveguide while in the asymmetric 
model some of the components of the continuum field are exponentially 
damped on one side of the waveguide. The physical interpretation of 
A (a;, v) will be discussed later. 

In the remaining interval, 7 4 = [ — k 2 K 2 , oo], the functions <o„ (n = 
1, 2, y) are all real and the functions M jl; {v) {j, k = 1, 2) have no 
poles. Therefore. 

dpM = p'M dv, (164) 
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where 

p' („) = i ( Wl + co 2 ) {col cos 2 («,«>) + w^ sin 2 (u u w)}/S>, (165) 

Pi»W = p'M = -«,(«i - toM + Wxtoa) sin («,.) COS («„«>)/©, (165) 

PmW = -«'(«i + w 2 ){«* sin 2 (u v w) + WiWj cos 2 («„w)}/SD, (167) 

©(„) = («£ + Wl o> 2 ) 2 sin 2 (2« y u>) + «.(«, + co 2 ) 2 cos 2 (2 W „w). (168) 

The spectrum for the TE fields of the asymmetric model consists of 
the (possibly empty) set of points v t , j = 1, 2, • • • , R and the interval 
— k 2 K x ^ v ^ oo. The transmitted field can now be written in the 
following way. 

«j°c».») = sit; iwwi"*iA^)fcfo)}«p i-iv^^i^,^ 

If''' exp {-iV^gjAC^^Sr^^Mp,^)/^)} dv 

7T J-k'Kt W-l J 

2 rtO 

J] / exp (-t'V-i'2)^, v)g k (y)p' ik {y) dv 

i,k-\ J-k'Kn 

+ E f exp { - VJfWft 9)*b)M) dv. (169) 

j, fc-l •'0 

The expression for e"'(x, 0) has been split up into a sum of parts in 
order to facilitate its physical interpretation. The first part represents 
the possible discrete, propagating modes which can be excited in the 
system. The form of these modes has been studied in detail elsewhere, 5 ,e 
and as pointed out earlier, unless condition (147) is satisfied, no such 
modes can be excited. In order to interpret the second term, consider 
within the waveguide in the region x < -wan incident plane wave 
of the form 

el 0) (x, z, v) = exp { -iV^vz - i^{v)x] . (170) 

At the surface x = —w, part of this wave will be reflected and part 
will be transmitted. However, at the surface x = w, the wave will 
suffer total internal reflection. The total electromagnetic field set up 
by e y 0) (x, z, v) is proportional to A(x, v) exp {-iV-vz}. The second 
term is then just a superposition of plane waves which are totally 
reflected at x = w. In Fig. 4 we give a schematic description of these 
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F J g- 4 T~^ schematic diagram of the totally reflected wave in the asymmetric 
step model. The wave is incident on the junction region at x = — w where it is 
partly reflected and partly transmitted. The partly transmitted portion is then 
totally reflected at x = w. 

waves. In microscopy the theory of the Becke line is based on just such 
a superposition of totally reflected plane waves. 10 The third term is a 
superposition of plane waves which are reflected and refracted at 
x = ±w. The last term is a superposition of waves bound to the surface 
2 = and propagating in the positive and negative x directions. 
The time averaged, transmitted power is 



P t = (2C0MO)" 1 Z V^ t | Vdpufyi) g l (v l ) 
1-1 

+ {hM/VTpM}gM | 2 

V^v | r.GOfcW + r 2 (u)g 2 (v) | 2 [p &)/&&)} dp 

k'K, 

+ (2^0)- f^ V~p 1 1) g^g^p'Mj du. (171) 

3.5 TM Fields For Asymmetric Step Model 

The TM fields for the asymmetric model present no new features, 
and we merely record the results. We have 

fo(x, v) = cos (K r a) x x), | x | ^ w (172) 

= cos (K r a) x w) cos {u 2 (x — w) } 

- faKt/usK,) sin {K r u> z w) sin {w 2 (x - w)) , x = w (173) 

= COS (KrU x W) COS {o} t (x + W)} 

+ (<a,KJ(a x K 9 ) sin (K r u x w) sin { Wl (a; + w)}, a; ^ -w (174) 
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fcfo v) = {KM) sin (K^x), | z | £ to (175) 

= (K V M sin (K r w z w) cos [ut(x - w) } 

+ (K a M cos (Xrto.w) sin {u 2 (x -w)}, x^w (176) 
= -(K B M sin (22> x u>) cos {o>i(a; + w)} 

+ (£,/«») cos (K T u x w) sin { Wl (s + «>)}, a ^ -w (177) 

where o>„(v) (n = x, 1, 2) are defined in (136) and A'„ and K f are de- 
fined in (105). Next, 

wi» = {(u x /K,) sin (JT^w) + »'(<*/**) cos (A^w)} 

■ { cos (K r o} x w) - i{u 2 KJo> z K 2 ) sin (K#,w) } -1 , (178) 
m-ceW = - {(«,//?„) sin (K,w x w) + i(«i/# x ) cos (X^,to)} 

• { cos (tf^w) - i{u l K t /u x K l ) sin (&>,«>) } _1 . (179) 
Then from (58) through (60) , (178), and (179) we obtain 

MM = N ik (v)/D(v) (j, fc = 1, 2), (180) 

where 

Nu(v) = -m - u^Kl/alKi^) 

+ (l + ^wtKj/wlKiXg) cos (2K rW;e Mj) 

- i(KJu z )MK x + w 2 //Q sin (21T rW ,w)], (181) 

AT 12 (v) = AT 21 (v) = (*/2)(«,/£, - «*/«, ( 182 ) 

- («J/Z3 + u^h/KtKt) cos (2K r o>,u>) 

+ i(ta x /K g )(f»i/Ki + «»/^») sin (2^^-w)], (183) 

D(„) = («,/£, + wtKJuXiK*) sin (2K r w,w) 

+ ifa/Kt + a> 2 /K 2 ) cos (2X r ft>,«;). (184) 

There are no points of the spectrum in Jj = [- oo, -k 2 K x ]. The 
only real poles of the functions M jk (v) are in the interval U_ = (-k 2 K x , 
-fc 2 A,). In Jo, <o x is real while on and wo are imaginary. If we let 

b(y) = KM*), P»(") = -*«»Wi n = *> 2 ' (185) 
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then the equation determining the poles, D(v) = is equivalent to the 
set of equations 

-v = k 2 K, + V \, -v = k 2 K 2 + v i , -v = k 2 K x - §= b 2 , 

(186) 
tan26w = (pJC./bKt + p 2 K,/bK a )/(l - p 1 p 2 K 2 c /b 2 K 1 K 2 ). 

In order that these equations have a solution, it is necessary and suf- 
ficient that - ° 

K r >K n (?i = l,2), (lg7) 

2wk{K,(K x -K 1 )/K,} i > tan -1 \KJvXK,- K 2 )/Kl{K x - K,))K 

If conditions (187) are satisfied, D(v) = has a finite number of real 
solutions m I 2 , vj, j = 1,2, • • • , S. 

If (186) and (187) are satisfied, we can write 

&r,»fo) = -N ik (vi)/D'(v t ), j, k = 1, 2, I = 1, 2, • • • , 5. (188) 
.lust as for the TE fields, it is true that 

{«*„(*)}' = IffM&rM* I = 1,2, •■■ ,S. (189) 

Because of (189) the functions ^(x, v t ) and \f/ 2 (.v, i> t ) appear in (49) 
for hl n (x, z) only in the combination 

*(X, v,) = V«<r, ,(»',) fcfo r,) 

+ (ftraW/VfanWI^C*, o,), i = 1, 2, • ■ • , flf. (190) 
If we define 

#o(3, v,) = VSff u (yj) cos (6(p,)x) 

+ [K&M/y/rM &("/)) sin (6(,,).t), (191) 
then because of (186) 

#(x, ",) = *»(*, pi), \x\ ^w (192) 

= ¥„(«;, *,) exp |p 2 (j/,)(w -a)}, 3 = «J (193) 

= * (-w, "/) exp [p^PjXw + x)\. x = -w (194) 
The remaining points in 7 2 are not in the spectrum. 

The remainder of the real axis, the interval —k 2 K x ^ v ^ °o forms 
the continuous spectrum. In the subinterval I 3 = \ — k 2 K, , — k 2 K 2 ), 
(ji x and oj, are real while o> 2 is imaginary. If we let 

K r <a x (v) = b(v), w,(i;) = p,(v), co 2 (v) = tjptOO, Kt/, , (195) 
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then we can write 

d*M = - IpM/KiWMvMv) dv (j, k = l, 2), (196) 

where 

8x(v) = cos bw + (p 2 K,/bK 2 ) sin bw, (197) 

8a ( v ) = -(p 2 /K 2 ) cos 6u> + (5/#.) sin bw, (198) 

AOO = {(&/JQsin26u> - (j? 2 /K 2 ) cos 26m;} 2 

+ {(pffltK./bKiK*) sin 26u> + (?,/#,) cos 2bw\ 2 . (199) 

When p e I« , ftfa r) and ^ 2 (x, r) appear in (51) for ftj°(*, z) only in 
the combination 

E(o>, *) = iiW^iC*, *) + «ito**C*, ")• ( 20 °) 

We have 

a(x, v) = cos { b(x — w) } 

- (p 2 K,/bK 2 ) sin {b(x - w)} , | x \ £ w (201) 
= exp {p 2 (w — x) }, x 2: w (202) 

= {cos 2bw + (p 2 K,/bK a ) sin 26u>} cos {p^x + w)\ 
+ (\/ Vl )\{bKJK z )sm2bw 

- (ptKj/Kt) cos 2bw\ sin {p t (x + w)}, x g> -to. (203) 
In the remaining interval, J 4 = [— MC 2 , oo], the functions w,, (n = 
1, 2, x) are all real and we can write 

cfcliW = ^00 cfr, (204) 

where 

(TfxW = - W*l + « a /X0 {(»'/«!) COS 2 (K^W) 

+ (fo l <o 2 /K 1 K 3 ) sin 2 (tf rWx w) } /», (205) 
■{■M = <rU») = " MK B )MK, - u 2 /K 2 ) 

■ [ol/K 2 , + u^/K^] sin (K^w) cos (K^w)/®, (206) 

riiW - I MK,)*(<*i/Ki + <»*/K 2 )[(fi»m/KiK%) cos 2 (K r a> x u>) 

+ ( Wl /A: e ) 2 sin 2 (X^,»)}/», (207) 



7T 
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JD = \MK g ) 2 + (w^/KJQ] 2 sin 2 (2K r u.w) 

+ MK^MK, + a> 2 /K 2 ) 2 cos 2 (2K& m w). (208) 

To summarize, the spectrum for the TM waves of the asymmetric 
model consists of the (possibly empty) set of points v t , I = 1, 2, • • • , S, 
and the interval —k a K l ^ v ^ oo. The transmitted field can be written 
as 

K'\x,z) =tt l^nWr'^^AWexp {-zV^ *)*(*,",) 
i-i *-i 

+ - ( " exp l-iV^z}Z(x, v) E*(")(P.W/^AW| rf" 

T J-k'K, ,=1 

2 »0 

+ 2 / exp {— «V^s}f/(s y i')AftG'Mft6')<b' 

l.t-l J-k'K, 
2 rtQO 

+ Z / exp (-V^ *}*,(*, ^*(")^*W^. (209) 

j.i-1 ■'0 

The time averaged, transmitted power is 

s 

P, = (2CO60)" 1 2 V-^,:| V<r u (v,) A x fr,) + \eMfVvtM\Wvi) ? 

V-y Is.W.W +s 2 {v)h 2 {v) | 2 {p^^/Zv.AW) <fr 

-Ifc'JCi 

+ (2coe ) _1 [ V^v £ h,{ V )*h k {vW ik {v) dp. (210) 



IV. APPROXIMATE SOLUTION OF THE INTEGRAL EQUATIONS 

In Section II we obtained general expressions for the reflected and 
transmitted fields for the TE fields in (18) and (27) and for the TM 
fields in (41) and (51). In (27) and (51) there appear the functions 
<Pi(x, v) and ypi(x, v) and the spectral density matrices p(v) and tr{v). 
A technique for determining these quantities in certain cases was 
illustrated in Section III by explicitly calculating them for the sym- 
metric and asymmetric step models. In order to complete the determina- 
tion of the reflected and transmitted fields, the functions &[ T) {1), 3C { u r) (l), 
g k {v), and h k {v) must be calculated. In Section II we showed that 
these functions were determined by the integral equations (30)-(31) 
and (54)-(55). 

We have been unable to solve these integral equations exactly for 
the general case. However, there are certain cases of great physical 
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interest, such as the electro-optic diode modulator, where excellent 
approximate solutions can be obtained. Let 

M n = max K n (x), ™» = min K n (x), (n = x, y, z) (211) 

and assume that 

(M n - m n )/m n « 1 (n = x, y, z). (212) 

Then the incident field impinges on an essentially uniform, plane 
dielectric interface, and the reflected field can be calculated as if the 
region z > were a uniform dielectric. Let R n (n = x, y, z) be suitably 
chosen, constant values for the dielectric tensor for z > 0. Then it is 
readily shown that for the TE fields 

&< r \l) = fi.(08i o (0, (213) 

and for the TM fields 

3C<°(0 = R ft (Z)3O0, (214) 

where the reflection coefficients are 

«.(0 = W) - M(WH»tf) + MKWP 1 . (215) 
R h (l) = {kMD - m/K)}{Km + nW.)}~\ (216) 

K = (R n ) h (n = x,y,z), (217) 

and a (I) is defined in (16). In this approximation, the total fields at 
z = for the TE and TM fields are, respectively, 

e v (x, 0) = ^ f^ r.(9s;°(Qe-"* dZ, (218) 

/>„(*, 0) = ^ /^ r fc (03ei o (i)e- <u di, (219) 

where the transmission coefficients are 

T n (l) = 1 + R n (l), n = e,h. (220) 

Now that e u (x, 0) and h u (x, 0) are known, g,(v) (j = 1, 2) can be cal- 
culated from (28) and h,(r) (j = 1, 2) can be calculated from (52), 
since e u (x, 0) = ej°(aj, 0) and fc,fo 0) = hl°(x, 0). 

We illustrate some features of the calculation of g k (i>) and h k (v) with 
the symmetric and asymmetric step models. We first note that if these 
models are used to study an electro-optic diode modulator, typical 
values of the parameters defining the dielectric tensors in (1) through (7) 
are 8 n = 3.31, A ^ 10~ 3 , 5 n ^ 2 X 10" 4 (n = x, y, z), A, = 0.96A, 
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A a = 1.04 A. Then M n - m n S 1.4 X 10" 2 , m„ « 10.9. Condition 
(212) is thus well satisfied. 

For the symmetric step model we let K n = K (n = x, y, z). If the 
functions S*°(0 and 3Cj°(Z) are sharply peaked about I = 0, then 
(21S) and (219) can be further approximated by 

e u ix, 0) = TM ~- r sj°(0r"- dl = r.(0)ei°(*, 0), (221) 

^,0) = r fc (0)*i°(a5, 0). (222) 

The calculation of g fc (v) and ft fc (v) is now reduced to quadratures. If 
the incident field is not sharply peaked, we define 

1 /** 
*A*i ") - ^ / */(*» ") c_ " x <**. (223) 

*,(*, ") = ^ J *,(*, v)[K{x) \ "V' ta dx, (224) 

so that 

0/00 = r r.(0Si°(0*y(/, ") ««, (225) 

*<fo - /" r»(03Ci°(0*,a, ") <*Z, J- 1,2. (226) 

If r is in the continuous spectrum, <*>;(£, v) and fy(Z, ?) are distributions 
which are easily determined with the aid of the relation 11 

/ e*" dx = 1 /(ia) + 7t5(<t), (227) 

where 8(a) is the delta function and when l/«r appears under an integral 
sign, it is assumed that the Cauchy principal value is taken. If v is in 
the point spectrum, *,-(/, v) and tf,(Z, v) are ordinary functions. 

For the asymmetric step model we letif n = \(K^ -f K 2 ), (n = x, y, z). 
For this model, a straightforward application of (28) and (52) fails 
in general if v is the point spectrum or if v e I 3 , because <pj(x, v) and 
fi(x, v) now grow exponentially as x tends to either plus infinity or 
minus infinity. This apparent difficulty is merely a reflection of the 
manner of convergence of the integrals defining g k (v) and h k (v). For 
our purposes here, it is enough to note from (169) and (209) that 
when v is in the point spectrum, the functions g k (v) and h k (v) do not 
appear independently, but only in the linear combinations 
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£ [spMr^piMgM = [ «.(*, o)*(x, *,-) dx, j = l, 2,---,R, 

t_l ■» — 00 

(228) 
E [8fuWr l M*W = f W*i 0)*(x, v,) dx, j = 1, 2,- • -, 5. 

1.= I ■»— CO 

(229) 

The integrals on the right of (218) and (219) are now well defined. 
Similarly, if v 1 1, , the relevant quantities to calculate are 

Z r k {v)g k {v) = f e„(x, 0)A(:r, r) dx, (230) 

E *W*tW = f UBi 0)E(aj, r) dx. (231) 

If v e I* , (28) and (52) can be applied directly. Now, all the techniques 
discussed in the case of the symmetric model can be applied here. 

V. SUMMARY 

In Section I we have defined a class of dielectric waveguide models. 
The waveguide is formed by an anisotropic, nonuniform dielectric 
filling the half space z > 0. The dielectric tensor is diagonal in the 
fixed coordinate system of Fig. 1, and the diagonal matrix elements 
are functions of x only, K n (x) (n = x, y, z). 

Integral representations for the incident, reflected, and transmitted 
fields were given in (15), (18), and (27), respectively, for the TE fields, 
and in (39), (41) and (51), respectively, for the TM fields. These rep- 
resentations are very general, holding for a large class of functions 
K n (x) and incident fields. These integral representations, however, con- 
tain the unknown functions <p,(x, v), xf/^x, v), p jk {v) and <r jk (v) (j, k = 1, 2), 
which are determined solely by the dielectric tensor, K n (x), and the 
unknown functions g k (v), h k (v), (k = 1, 2), 8^(0, and <JCj r) (0, which 
also depend on the incident field and the boundary conditions at z = 0. 
It was shown that this latter group of unknown functions are the solu- 
tions of two sets of integral equations, (30)-(31) for the TE fields 
and (54)-(55) for the TM fields. These equations are very complicated, 
and we have been unable to solve them exactly for any specific models 
of interest. 

In Section III we gave a detailed calculation of the functions <p t (x, v), 
fifa v), p ih (v), and a ik (v) (j, k = 1, 2) for both the symmetric and 
asymmetric step models. These calculations are important in their own 
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right, since the symmetric and asymmetric step models have been 
used extensively in the study of the electro-optic diode modulators. 1 " 6 
However, these computations also illustrate the technique for treating 
the whole class of piecewise constant models. This is important, for 
it is not yet completely established which is the correct model to use 
in exploring the behavior of the electro-optic diode modulator, and 
it is felt that any actual physical situation can be well approximated 
by a piecewise constant model. 

It should be noted that the success of the techniques used in this 
paper depends on being able to obtain exact analytic solutions of (21) 
and (43), or at least good analytic approximations to these solutions. 
There are a number of other models for which the exact solutions of (21) 
can be obtained, for example the continuous dielectric constant models 
described in Section III of NM. It is, however, much more difficult 
to find models, other than the piecewise constant models, for which 
(43) is solvable in terms of known functions. Nevertheless, the pos- 
sibility remains of investigating the TE fields for a fairly wide varity 
of models. 

The calculations of Section III provide a method of determining 
the discrete modes which is different from the methods used in earlier 
treatments. 6 ' 9 These calculations showed also that the asymmetry 
of the background light is accentuated in the asymmetric step model 
by total internal reflection at the junction region boundary. 

Finally, in Section IV it was shown that good approximations can 
be found for the functions g k (v), h k (i>), 3Cj r) (Z), and 8j r) (Z) in certain 
cases of physical interest. In particular, these approximations are valid 
for the electro-optic diode modulator. These approximations do not 
depend on a particular choice of the incident field. 

The final results of this paper then are integral representations for 
the fields for both the TE and TM fields. Of the various functions in 
the integrands, some have been determined exactly and good approxi- 
mations have been found for the remainder for a number of important 
models and for arbitrary incident fields. 

These integral representations are complicated in appearance, but 
when z is large enough, asymptotic expansions of them can be found 
which lend themselves to numerical analysis. In a subsequent paper 
asymptotic expansions of the transmitted fields will be presented for 
the symmetric and asymmetric step models in the case that the inci- 
dent field is Gaussian and numerical results for cases of experimental 
interest will be presented. 
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